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1 Introduction 

Perturbation theory is the standard method to study quantum field theories in 
the small coupling regime. However, the interplay of the perturbative expan- 
sion and the thermodynamic limit remains controversial. In particular, argu- 
ments were put forward that the infinite volume limit of perturbative coefficients 
does not give the correct infinite volume asymptotic perturbation expansion of 
asymptotically free theories In addition to the standard free (FBC), pe- 
riodic (PBC), and Dirichlet (DBC) boundary conditions for the spin model 
considered, a novel boundary condition was introduced in [P, namely superin- 
stanton boundary conditions (SIBC). The latter consist of Dirichlet conditions 
on the boundary of the system, and the additional freezing of one spin in the 
center of the sample. Perturbative coefficients were shown to have different 
thermodynamic limits for standard boundary conditions and SIBC. 

It was argued in ^ (see also [5]) that SIBC do not possess a well defined 
perturbation expansion, the third order coefficient being predicted to diverge in 
the infrared. Thus, perturbation theory was assumed to be consistent as the 
V ^ oo limit for standard boundary conditions is taken. 

This is a companion paper to ^ (so far the last contribution to the contro- 
versy, see citations therein), in which the volume dependence of perturbation 
theory coefficients for the 0{N) vector model with different boundary condi- 
tions was investigated up to third order, confirming the points of |2| regarding 
independence of the infinite volume perturbative coefficients for 'standard' b.c. 

The aim of this paper is to describe in detail the method used to compute the 
perturbative coefficients in and give a broader view of the results, including 
the IR divergence of SIBC correlators at third order. 
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2 O(A^) model: Perturbation theory 

The partition function for the O(A^) spin model on a 2-dimensional lattice A is 



Z : 



exp 



(1) 



where stands for a sum over nearest neighbour pairs. The delta functions 

constrain spins to have unit norm, i.e. S S S'^""'^ C Va; G A. 

The perturbative expansion of this model around the classical vacuum con- 



figuration Sx = •Si'^'' = (1, 0) Va; G A is constructed by writing 



Sx — [f^x, T^x), 



(2) 



and Taylor expanding in powers of tt . 

Assume for simplicity that the problem has no zero modes, i.e. the boundary 
conditions are such that the only vacuum configuration is {sl"''}. The case with 
zero modes will be treated in next section. 

The measure of the functional integral, in terms of tTx, is written as 



]ld^SxS{S'x-l) = l[^^=l^^ 



(3) 



with unconstrained tTx- 

The square roots can be exponentiated and added to the action terms in the 
exponent of This yields 



exp < ^ ^ (nx ■ ny + ^I-tt^ y/l-TT^) " \ 



Expanding in powers of /? ^ after a rescaling tTx — > tTx/^/P, 
Z^ |[nd--7r.]exp|-i ^,)^-^^(^^-,r: 



(4) 



vl 



(a;, v) 



16/3 



(5) 



(terms containing a single sum over x come from the change in the measure) 
2.1 Hasenfratz terms 

If the boundary conditions are such that there is a continuum of classical vac- 
uum configurations, obtained by rotation of S^'' (e.g. for FBC or PBC), the 
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corresponding zero modes have to be dealt with by introducing collective coor- 
dinates. As shown by Hasenfratz 5 , this amounts to adding an extra term to 
the action, 

-(7V-l)ln^Vl-7r2. (6) 

X 

The integrand of the partition function jSJl gets multiplied by 
f A^- 1 



exp • 



V a: X X y J 

(7) 

Notice that all terms in the exponent are suppressed in our notation by 
powers of V, the volume of the system. 

2.2 Vertices 

The quadratic part of gives a propagator S^^Gxy (superindices denote 'colour', 
subindices denote lattice points). The precise form of G^y depends on the 
boundary conditions, and will be discussed in next section. 
The vertices corresponding to lO and (0 are: 

• Order f3~^, two-?? vertex (o): 

V^:-U^-^)S''6.-^- (8) 



Order f3 ^, four-7r vertex (□): 

-xrkimrL ^ 

^ zwtu ^ 



(9) 



Order {3 ^, four-7r vertex (I 



zwtu 



2-^ 
V 

N- 1 



(10) 



• Order /3 ^, six-?? vertex (•): 



W^'Xt = ^ {-3n.<5,^,„p,(<5"<5""<5°^ + 14 terms) 



14 terms 



(11) 
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The meaning of the symbols used is the foUowing: is the number of 
nearest neighbours of site x. The value of Sx^x2...xn is 1 for xi = X2 = ■ ■ ■ — Xn, 
and zero otherwise. As for S^'j^', it is 1 if x and y are nearest neighbours, and 
zero otherwise (in particular, S^^- = 0). 

All terms explicitly dependent on the volume V are Hasenfratz terms, and 
should not be included in the perturbative calculations in absence of zero- 
modes. 



3 Propagator, boundary conditions 

The form of the 7f-propagator Gxy depends on the boundary conditions imposed. 
The quadratic form in the exponent of the integrand in (jSJ has the form 



with 



(12) 



(13) 



•^^x,y — ^x^xy ^xy i 

acting trivially on the internal O(A^) space. 

For DBC and SIBC, this quadratic form has no zero modes, and the matrix 
of propagators is just G = A4^^. However, for FBC and PBC, configurations 
with TTx = TToVx constitute the kernel of A4. Let P be the projector onto this 
space, and P-"- its orthogonal projector. Then 

M^MqP^ + OP, (14) 

where A^o is regular in the space of nonzero- modes. The matrix G is defined as 

G = Mo^P^+QP. (15) 

Let us discuss the form of M and G for the different boundary conditions 
mentioned. 



3.1 Free boundary conditions 

For definiteness, we work with a (strictly) 2-dimensional square lattice, with 
V = L X T sites, T = L. Rows are numbered from to T — 1, and columns 
from to L — 1. Lexicographically ordering sites in the lattice, and taking into 
account the different numbers of nearest neighbours lattice sites have, we write 
the V X V matrix ()13|l in terms of blocks of size L x L as 



/ a -1 
-1 b -1 
-1 b 



M 



-1 
a J 



(16) 
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where 1 is the L x L unit matrix, 



/ 2 



V 



-1 



-1 
2 



(17) 



and b = a + 1. 

Observe that ker is generated by (1, 1, . . . , 1)-^. 

Matrix Ai has a block tridiagonal structure, and its blocks are themselves 
tridiagonal matrices commuting with each other. By working within each eigenspace 
of the blocks, we can apply to the whole Ai an inversion procedure valid for 
tridiagonal matrices. 



If A is tridiagonal. 



A 



-1 



V 



fll 



-1 a2 



ar-2 

-1 



(18) 



its inverse is A -'^ = {^ab), with entries 

f-ab = AfcAb+i • • • Aa(l + AaAa+i(l + 
t^ab — l^ba^ 



■{l + \r-2K-i) ■■■)). a>b, 
a <b. 



(19) 



constructed from numbers Xj which can be computed recursively from the di- 
agonal elements of A: 



Ao — Oq 



Xj - [aj - Xj^iY 



1. 



(20) 



Possible divisions by zero can be avoided by rearranging the rows of A, as long 
as it is invertible. 

If our matrix Ai of eq. (|16|l were invertible, we could apply this construction 
to each of the eigenspaces of matrices a (eigenvalues a'^^\ eigenvectors v'^^\ 
h = 1, . . . , L — 1) and b — a + 1 (eigenvalues a*^''-' + 1, same eigenvectors). 
Procedure H20|l to compute the 'building blocks' A^'''' of the inverse matrix can 
be written, for eigenvalue a^^\ as 



A 



(h) 

T-l 



and the /^^^'^^ as 

l^ab - \ \+l ' 

(h) {h) 



= (aW + l-X)'^,) 



.Ar^(l + Ai''^Ai'i\(l + . 



T- 2, 



(21) 



■ (1 + Al'l^^Al'!.^) • • •)), a > b, 
a <b. 



(22) 
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Then, all we would have to do is construct G as 



L-l 



(23) 



h=0 



with P'^^^ = the L x L projector onto the h-th eigenspace of a, b. 

Now the eigenvalues of a (resp. b) lie in the interval [1, 3] (resp. [2, 4]), there 
existing just one eigenvector, w^"-* oc (1, 1, . . . , 1)^, with eigenvalue 1 (resp. 2). 
Then (|21|l can be carried out without problems if a > 1, but for the unit 
eigenvalue the last step is a division by zero, signalling the breakdown of the 
inversion procedure due to ^A being singular. 

We need to generalise the procedure to obtain G of the form (|15|l . To do 
this, first we add a regulator: 



i-VxV- 



(24) 



The eigenvectors of and b^ remain the same as for a and b, and their eigen- 
values are shifted by the small positive number e. Procedure H21|l can now be 
used for all eigenspaces. For eigenvalues strictly larger than 1 + e, the results 



are as before up to terms of order e. For the lowest eigenvalue a. 
Up, we obtain 



(0) 



Ai? = l-{j + l)s- 



(j + l)(j+2)(2j+3) 
6 



£2 + 0(e3), J = 0, 



X(0) _ ^ 



(T- i)(2r- 1) 



The elements of the inverse, /i^^^, , can be decomposed as 



,(0) 



1 



with 



(r-z)(r-i-») + T^-i 

2 2 6 ' 



1 + e of 
.,T-2, 

(25) 
(26) 

(27) 



Now it suffices to observe that the divergent piece is just the contribution 
coming from the zero mode (oc (1, 1, . . . , 1)"^) of as a whole. Indeed, it can 
be checked that the sums of the contributions of all other /i^'*-', and of Ji, to each 
row vanish. Then G is constructed as 



L-l 



which is the final result for FBC. 



(28) 



h=l 
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3.2 Periodic boundary conditions 

For PBC, fix — 4: at each site, hence 



/ c -1 
1 c 



M = 



-1 \ 



V -1 

in terms of blocks of size L x L, with 

/ 4 -1 

-1 4 



V -1 



c -1 
-1 c / 



-1 \ 



4 -1 
-1 4 / 



In this case, M can be immediately diagonalised by observing that 

T-l 



c — 2 cos 



with the T X T projector 



p(m) _ 1 m(a-&) 



LOT = exp z 



27r 



Now, of course. 



with the L X L projector 



L-1 

E 

n=0 



4 - 2 cos 



2Tm 

IT 



Q 



(n) 



n{c—d) 



271 

ujl = exp [ I — 



Inserting H33|l into H31() and using symmetry properties, we obtain 

T-l L-1 r 



4 - 2 cos 



771—0 n— 

with a total V x V projector 



27rr7i 



T 



2 cos 



27m 



L 



p(nin) 



COS [^m{a -b) + ^n{c - d)] 



V 



The contribution is zero, corresponding to the zero mode of Ai. This 
the decomposition p4|l . We can write directly G as in H15|l . 



G-EE' 



4 - 2 cos 



2TTm 
T 



2 cos 



27rn 

IT 



m— n— 

the prime meaning omission of the (m = 0, n = 0) term. 



p(mri) 



3.3 Dirichlet boundeiry conditions 



For 0-Dirichlet boundary conditions, spins along the boundary of the system 
are frozen to tTx = Vx € 9A. 

All Gxy with a: or y on the boundary thus vanish. This allows us to restrict 
the sums in the action and correlators (in terms of -k^) to the inner L x T = 
(L - 2) X (T - 2) lattice (matrices M and G will he V x V = {L - 2)(T - 2) x 
{L — 2)(T — 2)). All inner spins have 4 nearest neighbours, but the piece 
in A4 makes its structure different from the periodic case: 



/ d 
-1 



M = 



\ 



-1 
d J 



(38) 



in terms of blocks of size L x L, with 



/ 4 



V 



-1 

4 



4 

-1 



-1 
4 



(39) 



The construction applied above to FBC can be used in this case without the 
complication due to zero modes, since M. is regular here and all eigenvalues S^'^^ 
of d are 'safe'. The building blocks A^'*^ are computed as 



A, 



(h) 





(JC^))-!, 



(40) 



and the /x^^^ as 



f^ab — f^ba ' 



to get 



L-1 



h=0 



with P^^^ = and v^^^ the h-th eigenvector of d. 




3.4 Superinstanton boundary conditions 

Propagators for 0-superinstanton boundary conditions can be obtained from 
those for DBC by the following argument: 

Freezing of ttq = can be attained by adding a term ^ ttq to the action of 
the system with DBC, modifying the original quadratic form, 

M^Mx = M-XY, (43) 
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with Yab = SaoSbO: Sind taking the limit A —>■ +00 after aU calculations. 
Then G^^ = liniA^+oo Gx, with 



Gx = {M- XY)-^ = G(l - XYG)-^ = G + XGY 
In components, 



J2^"iYGYy 



.n=0 



(G\)xy — Gxy + XGxQ 

Hence, taking the limit, 



.n=0 



Gqj, — Gxy + 



XGxoGoy 
1 — AGoo 



YG. (44) 



(45) 



G^y — ^lini (GA)a:y — Gxy 



GxoG{ 



xO'-^Oy 



A — >oo 



G, 



00 



(46) 



4 The correlator 

We are interested in the 2-point function of Sx, whose perturbative expansion 
can be expressed in terms of tt as 



Sx ■ Sy) — 1 + ^ 
1 



1 



1 



1 
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We can write the perturbative expansion of {Sx ■ Sy) as follows: 

[s^-Sy)-i + ±f.+o{r^^^^^), 



with Ci = {Sx ■ S'y)^'-' polynomials in iV — 1, 



(47) 



(48) 



c, = ^c,J(7V-l)^ 1=1,2,... 

3 = 1 



(49) 



4.1 Feynman diagrams 

TTTT processes (Feynman diagrams) contributing to {Sx ■ Sy) to order l3~^, 
and are drawn in figure ^ 

The TTTTTTTr contribution comes from the diagrams in figure |21 (which should 
be drawn for each channel): 

Finally, figure |31 contains the unique Feynman diagram (up to channel re- 
ordering) representing the TrTTTrTTTrTr contribution to the correlator: 
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X — 2— Q — K 



K X 2 O K X O Q K 

Figure 1; Two-tt Feynman diagrams. 



Figure 2; Four-7r Feynman diagrams. 



4.2 Non-Hasenfratz terms 

Coefficients can be split into two terms, one coming from non-Hasenfratz 
contributions, the other term coming from Hasenfratz contributions: 



(50) 



Non-Hasenfratz contributions, in our notation, do not explicitly depend on the 
volume V. 

We now list all these contributions, up to and including order (i^'^ . Wo write 
Cij instead of c"j^' when the corresponding Hasenfratz contribution vanishes 
identically. 

We employ the following condensed notation: 



HI, = G. 



Gbc, 

^ab ~ GacGbd — GadGbc 



Jab — Gaa — Gbb, 

^ab ~ GacGbc ~ GadGbdi 



T)cd 



I /-<2 



(51) 



Order /3 ^ : there is only one contribution. 



Cll 2 {GxX "I" Gyy 2 Gxy ) 



o \ xy xy J 



(52) 



Order /3 ^: There are contributions proportional to iV — 1 and to {N — 1)^, 
C21 = — ^ -Ri^y — 2 -^^y ~'~ 2 {GzzH^y +GwwH^y —2GzwHxyHxy)j (53) 



Figure 3: Six-7r Feynman diagram. 
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and 



(54) 



Order /? there are contributions proportional to — 1, {N — 1)^ and 

(AT -1)3, 



C3I — 2 [^xx + ^xyi^xx + Gyy) 2 ^zwH^yH^y 

z, w 

— - y^(2 GzzH^J^ + GxxG^^ + GyyGy^ - 2GxyGxzGyz) 

— ^ I GxxGxzQVz ~ 2 GxyGxzQ^ + Gzw{Gzz + Gww)HxyH^y 

{z,w) 

+\ Kf - \ Pxf + (-3GL + GL)^fx\' + ^)] + (2: 2/)} 
+ 1 H {[imyHSyQi; + {GzpH^y-G^^H^yf + {z^w)]+{x^y)} 

{z, w),p 

~ ^ ^ ] 1^ (^-^xy^xyiGzzQ^p ~ GzwQ^p ~ GzpR^lj + G^qPz^) 
{z,w), {p,q) 

+ H^yGzp{H^xyQ'zl + H^yQll) + (p - q)) +{z^ w)] + {x ^ y)]. 
The 0{(}~^) contribution proportional to {N — Vf is 

^32^' = - g [^(G^x + Gyy) - {GxxGyy + 2Gly){Gxx + Gyy) 



(55) 



2^12 

xy 



- g X/ { [ ^ 2 ^^"^ ~ ^"''^ ~ ^ ~ ^ GzzGww)H^y + ^ -R: 

(2, to) 

+4 GxxGj;2(G2,2 — Gj^^) — Gzw{Gzz + Gww)HxyHxy 

-2 G,„G,,Q^^ - 2 G,,G,^Q^|i + (^ ^ w)] + {x ^ y)] 
+ ^ E { [2 HiyHlyGzpJz^ + (G^^ - G^p)//|,2 + (^ ^ w)] + 



2/)} 



{z,w) ,p 



pq 
zw 



{zyw), {p, q) 



~l~ ^xy^xyi^ zp(.-^^w ~^ JzwJpq H~ ^zz^pp ^wwGqq) 



J zw ^pq ^zq Jpq ^zw ^wp] 



-\-{p ^ q)^ -\- {z ^ w) -\- {x y)Y 



(56) 
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Finally, the 0{(3 ^) contribution proportional to {N — 1)^ reads 
1 



„n.H. _ 



~ 1C "^Xvi^XX + Gyy) 



{z,vj) 



+{p ^ qij + {z ^ w) +{x<r-^y)Y 



{z,w), (p,q) 



(57) 



4.3 Hasenfratz terms 



Hasenfratz terms, in our notation, come suppressed by powers of the volume V 
of the system. They should be included for PBC and FBC, to deal properly 
with the rotation zero-mode of the vacimm manifold. 

We list all Hasenfratz contributions which do not identically vanish. 

Order there is a contribution proportional to (A'' — 1)^. 



xy 



(58) 



Order j3 contributions to C32 and C33. 
C32 = ^ {pzzH^y + GxxG'^^ + GyyGy^ — 2 GxyGxzGyz^ 

z 

+ ( y ) X/ GzwHxyH^y 

^ ^ z,w 

~ 2V {^-^xyi^xyQzp ~ ^xyQwp) (GzpHxy — GwpH^y) |, 



(59) 



and 



^fi — \pzzHxy + Jxy{Gxz ^yz)} 

~^Qj/^ ^G ZzH xy ~\~ GyjyjH^y — AG ZW H xyH ^y^ 



G^p){H^y H^y )+ 2 JzwHxy{GzpH^y GyjpH^y)^. 

(60) 



(2, If), p 
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5 Results 



C programs were written to compute the coefficients given by the above expres- 
sions for square lattices of sizes 2x2 through 120 x 120. 

The main results were reported in 4 , to wit: standard boundary conditions 
give rise to coefficients agreeing with each other in the infinite volume limit, 
while SIBC coefficients disagree with them in this limit, actually diverging at 
third order, as predicted in [2]. 



5.1 Standard boundary conditions 

Figures01through|51compare coefficients ci, C21, C22, C31, C32 and C33 for standard 
boundary conditions: FBC, PBC, DBC. Data is plotted only for lattices larger 
than 20 x 20 for ease of inspection, and the region L > 100 is showed separately. 
The perturbative expansions for these boundary conditions can be shown to 
agree coefficient by coefficient in the thermodynamic limit, with limit values 
agreeing with |S]: 



F,P,D L- 



F,P,D L- 
"21 



1 

4' 
1 

32' 



F,P,D L-^cx „ 
C22 *■ U, 



F,P,D L^oo _ riri707 
C^l > —0.00727, 



F.P,D L^oo „ 

C32 — > -0.006, 

F.P,D L-»oo „ 

C33 — * 0- 



(61) 



5.2 SIBC: different thermodynamic limit 

SIBC coefficients up to second order are plotted in figurcEH Their convergence 
is extremely slow compared with the previous cases. Yet their behaviour seems 
compatible with the conditions 



„si 
'-21 



„SI 



L — *oo 



1 

4' 
1 

32' 



(62) 



necessary for the agreement, as L ^ 00, of the perturbative expansion of the 
Abelian {N = 2) model with the corresponding expansion for standard boundary 
conditions. This agrees with the results of 1 , which used another method. 



5.3 SIBC: infrared divergence at third order 

Third order coefficients for SIBC are plotted, this time as a function of logL, 
in figure ^2 Both C32 and C33 are seen to diverge (logarithmically in L), a phe- 
nomenon predicted in [2]. The behaviour of C31 is unclear, but the curve also 
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-0.249 



-0.2495 




125 



Figure 4: ci as a function of L for FBC, PBC and DBC, and blow-up of the L > 100 
region. 
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Figure 5: C2i as a function of L for FBC, PBC and DBC, and blow-up of the L > 100 
region. 
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Figure 6: C22 as a function of L for FBC, PBC and DBC, and blow-up of the L > 100 
region. 
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Figure 7: csi as a function of L for FBC, PBC and DBC, and blow-up of the L > 100 
region. 



17 




18 



0.0003 



0.00025 



0.0002 



0.00015 



0.0001 



5e-05 



-5e-05 



2.5e-05 



2e-05 



1 .5e-05 



1e-05 



5e-06 



-5e-06 



-1e-05 




125 



Figure 9: C33 as a function of L for FBC, PBC and DBC, and blow-up of the L > 100 
region. 
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Figure 10: ci, C2i and C22 as a function of L for SIBC. 
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Figure 11: C31, C32 and C33 as a function of logL for SIBC. 
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suggests an infrared divergence. Thus, the (3~^ coefficient presumably diverges 
logarithmically in the infrared for all N > 2. However, the Abelian case N = 2 
could yet prove convergent if the divergent parts of C31 , C32 and C33 cancel — 
in this case it would be interesting to check whether they agree with the corre- 
sponding coefficient with standard boundary conditions in the thermodynamic 
limit. 

6 Conclusions and outlook 

The perturbative expansion of observable {Sx • Sy) for the vector O(A^) model 
has been computed up to order for different boundary conditions. Consis- 
tency of the usual infinite volume limit of the perturbative coefficients has been 
checked up to this order for standard boundary conditions (PBC, FBC, DBC). 
Divergence of the third order SIBC perturbative coefficient in the infrared limit, 
as predicted in ^ , has been explicitly shown. It would be interesting to extend 
the analysis to larger lattice sizes (the current upper limit L = 120 was dic- 
tated by array storage requirements in available computers): the behaviour of 
the total coefficient could have a thermodynamic limit for the Abelian case 
N = 2. 

Expressions for Feynman diagrams have been computed exactly, and the tt 
propagator has been obtained exactly for all boundary conditions considered, 
using a method to diagonalise block tridiagonal matrices — computations being 
carried out by means of C programs. While the usefulness of this method 
depends critically on the particular boundary conditions used, a generalisation 
for other kinds of boundary conditions being thus unlikely, it is particularly well 
suited for exact computations in 2d rectangular systems. It could for instance 
prove useful to analyse the expansion of these models. 

The long standing problem of the correct perturbative expansion of asymp- 
totically free quantum field theories in the thermodynamic remains, all in all, 
unsolved. However much circumstancial evidence is gathered for the ortho- 
dox approach, for instance in the course of this work, a proof of the rule that 
standard boundary conditions provide the correct asymptotic series for these 
theories would be needed (which should involve a nonperturbative definition of 
the theories, since 'asymptotic' series need a definite function to be asymptotic 
to!). 
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